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Abstract. We examine the near collapse dynamics of a self-gravitating 
magnetized electron gas at finite temperature, taken as the source of a Bianchi- 
I spacetime described by the Kasner metric. The set of Einstein-Maxwell 
r—t field equations reduces to a complete and self— consistent system of non-linear 

autonomous ODE's. By considering a representative set of initial conditions, the 
numerical solutions of this system show the gas collapsing into both, isotropic 
("point-like") and anisotropic ("cigar-like") singularities, depending on the 
intensity of the magnetic field. We also examined the behavior during the 
collapse stage of all relevant state and kinematic variables: the temperature, the 
expansion scalar, the magnetic field, the magnetization and energy density. We 
notice a significant qualitative difference in the behavior of the gas for a range of 
temperatures between the values T ~ 10 3 K and T ~ 10 7 K. 



PACS numbers: 98.80.-k, 04.20.-q, 04.40.+-b, 05.30.Fk, 95.36.+X, 95.35.+d 
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1. Introduction 

Astrophysical systems provide an ideal scenario to examine the effects of strong 
magnetic fields associated with self-gravitating sources under critical conditions. 
In such conditions, we expect non-trivial coupling between gravitation and other 
fundamental interactions (strong, weak and electromagnetic), and from this interplay 
important clues of their unification could emerge. 

The presence and effects of strong magnetic fields in compact objects (neutron, 
hybrids and quark stars) have been studied in the literature (see [TJ [21 El HI El El El [S] , 
and references quoted therein), assuming various types of equations of state (EOS) 
has been obtained and considering in some of these papers numerical solutions of the 
equilibrium Tolman-Oppenheimer-Volkov (TOV) equation. 

As proven in previous work [TJ [H El EH HI], the presence of a magnetic 
field is incompatible with spherical symmetry and necessarily introduces anisotropic 
pressures. However, by assuming that anisotropics and deviations from spherical 
symmetry remain small, several authors [SJ 1101 112] have managed to compute 
observable quantities of idealized static and spherical compact objects under the 
presence of strong magnetic fields by means of TOV equations that incorporate these 
anisotropic pressures. Moreover, it is evident that much less idealized models would 
result by considering the magnetic field and its associated pressure anisotropies in the 
context of TOV equations under axial symmetric (or at least cylindrically symmetric) 
geometries [5]. 

As an alternative (though still idealized) approach, and bearing in mind the 
relation between magnetic fields and pressure anisotropy, we have examined the 
dynamics of magnetized self-gravitating Fermi gases as sources of a Bianchi I space- 
time [13l [HI [15], as this is the simplest non-stationary geometry that is fully 
compatible with a the anisotropy produced by a magnetic field source. 

Evidently, a Bianchi I model is a completely inadequate metric for any sort of a 
compact object, as all geometric and physical variables depend only on time (and thus 
it cannot incorporate any coupling of gravity with spatial gradients of these variables). 
However, the use of this idealized geometry could still be useful to examine qualitative 
features of the local behavior of the magnetized gas under special and approximated 
conditions. Specifically, we aim at providing qualitative results that could yield a 
better understanding of the conditions approximately prevailing near the center and 
the rotation axis of less idealized configurations, where the angular momentum of the 
vorticity and the spatial gradients of the 4-acceleration and other key variables play 
a minor dynamical role. 

The main objective of the present paper is to include the effect of the temperature 
in the study that we accomplished in [13.. We aim at addressing the question 
of whether a finite temperature produces a significant dynamical effect in "slowing 
down" or reversing the evolution of the magnetized gas in the collapsing regime, 
as such qualitative difference may be related to the stability of the self-gravitating 
configuration. We also analyze the relations between the magnetization and magnetic 
field, and the energy density and temperature. 

The paper is organized as follow. In section II we derive the EOS for a dense 
magnetized electron gas at finite temperature. In section III we lay out the dynamical 
equations for the evolution of our model by writing up the Einstein-Maxwell system 
of equation for the specific source under consideration. Einstein-Maxwell equations 
are written in section IV as a system of non-linear autonomous differential equations, 
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rewriting it in section V in terms of physically motivated dimensionless variables. The 
numerical analysis of the collapsing solutions and the discussion of the physical results 
are given in the section VI. Our conclusions are presented in sections VII. 

2. Magnetized Fermi gas as a source of a Bianchi I background geometry 

Homogeneous but anisotropic Bianchi I models are described by the Kasner metric 



ds 2 



-dt 2 + Q x {if dx 2 + Q 2 (t) 2 dy 2 + Q 3 (t) 2 dz 2 , 



(1) 



so that spatial curvature vanishes and all quantities depend only on time. Assuming a 
comoving frame with coordinates x a — [t,x,y,z] and 4- velocity u a — 5f, the energy- 
momentum tensor for a self-gravitating magnetized gas of free electrons is given by: 

T b a = (U + P)u a u b + P8 b a + Il b a , P=p -'^ i ( 2 ) 

where B is the magnetic field (pointing in the z direction), U is the energy density 
(including the rest energy of the electrons), M is the magnetization of the gas, P is 
the isotropic pressure and Ii b is the traceless anisotropic pressure tensor: 

n fc a = dm 5 [n,n,-2n,o], n = -^-, (3) 

We can write the energy-momentum ^ tensor as: 

T, Q = diag[-[/,P ± ,P ± ,Pj|], (4) 

which respectively identifies P± and Py as the pressure components perpendicular and 
parallel to the magnetic field. Notice that the anisotropy in T b a is produced by the 
magnetic field B. If this field vanishes, the energy-momentum tensor reduces to that 
of a perfect fluid with isotropic pressure (an ideal gas of electrons complying with 
Fermi-Dirac statistics) . 

The equations of state for this magnetized electron gas can be given in the 
following form [16] : 



Pi 



-n = \f3 

A Th {P,n,4>) > 

OO s i 

P ± = X(3 2 J2nC 1 ' ^ 
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The particle number density is 

n = —fi , 

771 e 2 
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-C 4 (<M + X; 
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(f3,/j,,4>) 



with a n = vT-j- 2n/3 and Ci, C 2 , C3, C4 are given by: 

00 

f 1 dx 



(5) 
(6) 

(7) 
(8) 



Gravity induced evolution of a magnetized fermion gas with finite temperature 4 



,.2 



fix 



C 2 (0, M ) = / , — - -j^=, (10) 
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C 3 (<M = I y -^-^^dx, (11) 
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C 4 (0,m)= / N <*c- (12) 



o 



l + exp(^±&) 

The equations ((SJ and ([5]) can be accommodated as [HI HZ] : 

Pl = P|| - BM, (13) 
where M, the magnetization, can be written in the following form: 

M = M Y M {[3, l iA), (14) 
with Mq = 2tthb I (A 3 .) anc ^ -^M * s obtained of the expressions (JSJ) , ([5]) and ([1 



G^i/ — R^v "^RQfiv — kT^ u , (15) 



In the previous expressions \x is the dimensionless chemical potential normalized 
by the rest energy, 4> — kT/m e , j3 = B/B c , where B c = 4.414 x 10 13 G is the critical 
magnetic field, is the Bohr magneton and A = 87rm e /(A 3 ) where A c is the Compton 
wavelength of the electron. 

In Appendix A the integrals (IH))- (IT2l have been transformed into equivalent 
expressions in order to facilitate the numerical calculations. 

3. Einstein Maxwell equations 

Since we are interested in the critical relativistic regimen, the dynamics of the 
magnetized gas whose EOS we have described in the previous section must be studied 
through the Einstein field equations in the framework of General Relativity: 

1 

2' 

together with the balance equations of the energy-momentum tensor and Maxwell's 
equations, 

T^ ;v = 0, (16) 
, v = 0, F [fW . a] = 0, (17) 

where k = &t:Gn and Gat is Newton's gravitational constant, while square brackets 
denote anti-symmetrization in (|17[) . 

Assuming absence of annihilation/creation processes, so that particle numbers 
are conserved, leads to the following conservation equation: 

n% = 0, n a = r 1U a , (18) 

where rj es the particle number density. From the field equations (|T5|) we obtain: 

_qx ^ Q2Q3 _|_ Q2 _|_ Q3 _ _ kP± / lg \ 
Q2Q3 Q2 Q3 

y Q1Q3 Q\ Q3 
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-G/ = g# + £ + £ = -^|, (2D 

_^ f _ Q1Q2 + Q1Q3 + Q2Q3 _ ^2) 
Q1Q2 Q1Q3 Q2Q3 

where Q = Q ]a u a = Q,t- From the conservation of the energy-momentum tensor (|16p 
we obtain: 

*-(| + i)^ + ^-|rt +l o- < 23 > 

Maxwell's equations (fT7|) yield: 

and from the particle number conservation (|18[) leads to: 



* (25) 

4. Local kinematic variables 

Einstein-Maxwell held equations are second order system of ordinary differential 
equations (ODE's). In order to work with a hrst order system of ODE's, it is useful 
and convenient to rewrite these equations in terms of covariant kinematic variables 
that convey the geometric effects on the kinematics of local fluid elements through 
the covariant derivatives of u a . For a Kasner metric in the comoving frame endowed 
with a normal geodesic 4-velocity, the only non-vanishing kinematic parameters are 
the expansion scalar, 6, and the shear tensor a a p: 

© = U a - a , G a p = M( 0;( 8) — — h a p , (26) 

where h a /3 = u a up + g a p is the projection tensor and rounded brackets denote 
symmetrization. These parameters take the form: 

Wl W2 H3 

af = diag [a*, 0* , af, 0] = diag [£ x , £ 2 , S 3 , 0] , (28) 

b^c (a,b,c= 1,2,3). (29) 



where: 



3 Q a 3 Qb 3 Q c ' 

The geometric interpretation of these parameters is straightforward: represents the 
isotropic rate of change of the 3-volume of a fluid element, while erg" describes its rate 
of local deformation along different spatial directions given by its eigenvectors. Since 
the shear tensor is traceless: c Q Q = 0, it is always possible to eliminate any one of the 
three quantities (E 1; £ 2 , £3) in terms of the other two. We choose to eliminate £1 as 
a function of (£2, £3)- By using equations (|2"T)l and (f2"9"|) we can re-write the second 
derivatives of the metric functions in (fT9"|) . (f2"U|) y (|2ip as first order derivatives of 0, 
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£ 2 and £ 3 . After some algebraic manipulations it is possible to transform equations 
P9P ~ P5|) as a first order system of autonomous ODE's: 

S 2 = |(P x -P||)-eS 2 , (30a) 

B = 2B (x 3 - , (30c) 

A = 77-7 [/i*?,r - / 2 t/ r ] , (30d) 
Det 

f 

Det 

together with the following constraint: 



^ = — W,p - hvj , (30c) 



^[/ = -(£ 2 ) 2 -£ 2 £ 3 -(£ 3 ) 2 + ^, (31) 



where: 



Det = U^.t - U lT r) <fl , (32) 
/1 = (Sa - ^) (Pi + CO - (e 3 + I) (P\\ +U)~ U <B B, (33) 

h =Qv + V,bB. (34) 

These first order equations form a complete and self-consistent system whose numeric 
integration fully determines E 2 , £3, B, fi, T, and thus allows us to study the dynamical 
evolution of a local volume element of a gas of magnetized electrons. 

5. Dynamical equations 

By introducing the following dimensionless evolution parameter, 

H - 6 d - 1 d (35) 

3 ' dr Hq dt ' 

together with the dimensionless variables, 

H=^r, S 2 = -^-, S 3 = ^-, P = 1T> A* = — , (36) 
n tl no t5 e m e 

where we have denoted by fl the usual chemical potential and Hq is a constant with 
inverse length units that sets the characteristic length scale of the system, which we 
have chosen as 3H% = k\ H = 0.86 x 10~ 10 m -1 , so that l/H = 1.15 x 10 10 m is 
of the order of magnitude of an astronomic unit. It indicates that our simplified model 
is examined on local scales smaller than cosmic scales. In cosmological sources and 
models [T5] Hq — 0.59 x 10 _26 m _1 would play the role of the Hubble scale constant, 
this value is a much greater length scale. The functions 5* 2 and S3 are the components 
of the shear tensor normalized with this scale, while r is the dimensionless time. 
Substituting (|3"f?|) into the system (|30aj) - (|30e[) we obtain: 

5 2 , r = Dl - r,| - mS 2 , (37a) 

5 3 , T = 2 (r,| - r ± ) - 3H5 3 , (376) 
0, T = 2(3 (S 3 - 2H) , (37c) 
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Det 

while the constraint (1311) becomes 



r = ^(r^/i + r^/ 2 ), (37 d) 

^(rj + rjj). (37c) 



ST V = -Si -Si- S 2 S 3 + 3H 2 , (38) 
and the auxiliary parameters of the previous system take the form: 

fx = (S 3 - 2H) (T ± - 2T UtP P) - (S 3 + U) T\\ - 3TuH, (39a) 

h = mr v + 2r„,0 (s 3 - m) , (396) 

Det = Tu,nT Vt 4, - Tu^T^. (39c) 

In the following section we undertake the numerical study of this system, focusing 
concretely in the the collapsing regime. 

6. Numeric analysis and physical interpretation 

Since the energy-momentum tensor (jlj) takes the perfect fluid form (P|| = Pj_) for zero 
magnetic field, we can identify the magnetic field as the factor introducing anisotropy 
in the dynamical behavior of the fermionic gas. In particular, this anisotropy in the 
stress (Pn 7^ P±) must yield different evolution in different directions, which must 
be evident in a critical stage such as the collapsing regime. Intuitively, we expect 
an isotropic point-like singularity if the pressure is isotropic, as pressure diverges in 
all direction, but a large pressure anisotropy (which necessarily corresponds to large 
magnetic field) should lead to a qualitatively different direction dependent critical 
behavior of the pressure that should result in an anisotropic cigar-like singularity 
characterized by the divergence of only the pressure parallel to the magnetic field. 
As we show in figure [TJ the pressure parallel and perpendicular for different initial 
temperatures does exhibit the expected behavior: for initial conditions of small 

anisotropic evolution Pn — > oo with P± — > that signals a cigar-like singularity (the 
case Pi I = Pj_ with zero magnetic field is shown as a comparison). 

Since V = ^f— det g a p = Q1Q2Q3, we obtain by means of pTj) and (|35|) the local 
volume in terms of H : 

V (t) = V (0) exp ^3 J UdT j , (40) 

where we remark that the sign of % (r) implies expansion if % (r) > 0, and collapse if 
T-L (t) < 0. Besides this point, equations (|27]l and ([29)) lead to: 



temperatures 0(0) = 10 7 (~ 10 3 K) we have Pm w Pj_ and a point singularity: 
Pi, P_l — > 00, but for a larger initial value 0(0) = 10 -4 (~ 10 6 K) we have the highly 



i(t) = Qi(0)exp 



[U + Si)dT 



T = 



(i = 1,2,3). (41) 



where S\ = — (S% + 53). 

For the numerical study we assume a magnetized electron gas at high density: 
A*e(0) = 2, which means that the chemical potential is 2m e . The initial values of the 
magnetic field and temperature were chosen in the ranges /3(0) ~ 10~ 5 to /3(0) ~ 10 -4 
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Figure 1. Pressure parallel and perpendicular to the magnetic field. We plot 
P||/A vs P±/X for initial temperatures 0(0) = 10 — 7 , 10 — 4 where A = 87rm e /(A 3 ). 

The case j3 = is shown for reference. 

and (f> ~ 10~ 7 to ~ 10~ 3 respectively Q Together with "H (0) < 0, we consider the 
conditions S^O) = 0, 53(0) = 0, +1, which correspond to the cases with zero initial 
deformation and initial deformation (shear) in the direction of z axis respectively. The 
calculation has been done using the fourth-order Runge-Kutta method with the local 
truncation relative error less than 10~ 6 . 

The numerical solutions for the function T~L for the assumed values H(0) show that 
% — > — 00, which implies that the volume element evolves to a singularity (see equation 
l|40p). This is exemplified in figure [3J where numerical solutions are displayed for the 
expansion scalar. These curves correspond to different values of the initial temperature 
in the range cf>(0) — 10~ 7 to (f)(0) — 10 -3 , fixing the rest of the initial conditions on the 
values n(0) = 2, /3(0) = 5 x 1CT 5 , S 2 (0) = and 5 3 (0) = 1. Notice that in this regime 
(as given by these initial conditions) where the high densities are dominant, we do not 
obtain a direct relation between the values of initial temperatures and the collapse 
time. In all the configurations the strength of the magnetic field diverges during the 
collapse of the volume element: this happens independently of the initial conditions. 
However, the collapsing time diminishes when the values of the initial magnetic field 
increases, which agrees with the results obtained in [13] . where T = was assumed. 

The expression (j43~|) shows how the evolution of the terms Si + % to ±00 implies 
that the metric coefficient Qi evolves to either +00 or 0. This evolution is characteristic 
of an anisotropic "cigar-like" singularity, since two metric coefficients evolve to zero 
and the third one to +00, whereas when all metric coefficients tend to zero, the 
singularity is isotropic "point -like" . See the definitions of these types of singularities 

\ Let us remark that <f> = kT/m e , f3 = B/B c , with B c = 4.414 X 10 13 G and fi is the chemical 
potential normalized /i e = /j, e /m e , the magnetic field is in the range B ~ 10 8 G to B ~ 10 9 G and the 
temperature varies between T - 10 3 K and T ~ 10 7 K. 
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in [H]. 

For initial conditions with zero shear: 5i(0) = 5*2(0) = 5*3(0) = 0, we 
always obtain a point-like singularity, independently of the selected values for the 
remaining initial conditions. However, when the initial deformation is positive in 
the direction of the magnetic held (z axis), we can obtain either cigar-like or point- 
like singularities. The nine curves displayed in figure [3] correspond to three sets of 
functions: Si + H (i = 1,2,3), each one corresponding to these different values 
of the temperature: (f)(0) = 1(T 7 , 10~ 6 , 10~ 3 . All curves were obtained taking as 
initial conditions: fi(0) = 2, /3(0) = 5 x 10~ 5 , 5 2 (0) = and 5 3 (0) = 1. For 
the higher initial values of temperature (f> (0) = 10 -6 , 10~ 3 the gas collapses into a 
cigar-like singularity in the direction of the magnetic field, whereas for the lower value 
4> (0) = 10~ 7 we obtain a point-like singularity. As expected, an anisotropic singularity 
occurs as the magnetic field increases, even for relatively low initial temperatures if 
the initial magnetic field (3(0) is sufficiently large. 

For initial values in the range 0(0) - 10~ 7 to 4>(0) ~ 10~ 4 the temperature 
decreases as the collapse proceeds, although grows up quickly for values near the 
collapse time, while for (f>(0) ~ 10~ 3 and higher values it increases. 

We display in figure 0] two dashed dashed curves for the evolution of the 
temperature from the initial value (f> (0) = 10~ 3 , together with solid curves that depict 
temperatures starting at (f> (0) = 10~ 7 . In both cases the initial conditions correspond 
to zero initial shear and positive deformation in the z axis. The initial values of 
magnetic field and chemical potential were fixed at (3(0) = 5 x 10~ 5 and fi(0) = 2 
respectively. Notice that the symmetric configuration delays the collapse (yields a 
longer collapsing time), which may indicate a connection with stability of local fluid 
elements in compact objects. This could provide an important clue on the stability of 
a compact object made of a dense magnetized gas. However, verifying this possibility 




Figure 3. The plots of the functions Si + "H for i = 1,2, 3. The system has an 
initial deformation (shear) in the same direction of the magnetic field along the 
z axis. The collapse is in the form of a "cigar— like" singularity in the z direction 
for 4>{0) = 10~ 3 , 1CT 6 and a "point-like" singularity for <f>(0) = 10 -7 . 



is beyond the scope of this article. 

We examine now the behavior of the magnetization and the energy density of the 
system during the collapse. We remark that the r-depending functions: ft, /i, <j) follow 
from the numerical solution of the system (I37a[) - (l37ep and (|38p . which allows us to 
compute the magnetization M (f3, /i, <fi) and U ((3, /i, <fi) for all values of r. We have 
depicted curves of M vs. the magnetic field and the energy density vs. temperature 
H for initial conditions: /3 (0) = 5 X 1CT 5 , S 2 (0) = S 3 (0) = 0, and for two different 
values of the initial temperature cf> (0) = 10 -7 , 10~ 3 . These curves are shown in figures 
Eland© 

Notice how at low temperature values the magnetization has an oscillatory 
behavior with respect to the magnetic field, which correspond to the well known 
Haas- van Alphen effect [20]. In contrast, for high initial temperatures, the system is 
heated as the collapse proceeds, with the magnetization behaving monotonically with 
respect to the magnetic field. These results agree with the study performed in |21) . 
where the authors expressed the magnetization of an electron gas at finite temperature 
as the sum of two terms: one oscillatory and the other monotonic, which predominate 
at respectively low and high temperatures. 

On the other hand, the energy density increases with the temperature for the 
initial value 4> (0) = 10~ 3 , which does not occurs for the initial value (j> (0) = 10~ 7 . 

The behavior of the magnetization and the energy density can be explained by 
taking into account that for the initial value 0(0) = 10~ 7 , as we discussed above, 
the system is cooled (temperature decreases) during the collapse process with an 
increasing magnetic field. This is an indication that the effects of the magnetic field 



§ Remember that M and U can be written in the following form: M = MqTm and U = Ar^. 
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Figure 4. Temperature of the system. The left hand side of the graph displays 
the temperature values (0i) corresponding to the two solid curves. The right 
hand side displays the temperature values (02) corresponding to the two dashed 
curves. 




Figure 5. Magnetization versus magnetic field. The axis of the right hand side 
(Tm2) depicts the magnetization (the solid curve), while the axis of the left hand 
side (rjifi) corresponds to the dashed curve. 
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Figure 6. Energy density versus temperature. The temperature values of the 
curve draws with solid line should be observed in the lower axis (<p2) and in the 
top axis (<pi ) , those of the curve drawn with dashed line. 



(alignment of electron spins with the magnetic field) predominate over the effects of 
temperature (random motion) in the system, a situation in which the Haas-van Alphcn 
oscillations appear in the magnetization curves, with energy density decreasing with 
temperature in the collapse process. On the other hand, for high initial temperatures 
(<fi (0) = 10~ 3 ), the system is heated (temperature increases) during the collapse, 
which shows the predominance of the effect of temperature over the effect of the field, 
and showing the magnetization and the energy density following a monotonic behavior 
with respect to the magnetic field and the temperature, respectively. That is to say, 
when (j) (0) ~ the temperature of the system decreases in the collapse process, 

while the magnetic field increases and the gas evolves to strong magnetic field regimen, 
where the effect of the temperature is negligible in a first approximation. This is not 
the case for </>(0) ~ 10~ 3 , for which the temperature rises significantly during the 
collapse and its effects are significant. 

It is worthwhile commenting on the relation between the temperature, magnetic 
field and magnetization in magnetized gases. In an "earth bound" context when 
the gas is not self-gravitating, the magnetic field acts as an "external" agent whose 
effect is to increase the magnetization by aligning the magnetic moments of the 
electrons. For higher temperatures the internal energy increases and there is more 
resistance to this effect, resulting in a relatively low magnetization with the electron 
magnetic moments alignment being relatively random, while at lower temperatures 
the opposite effect occurs: magnetization is high and magnetic moments strongly 
align with the magnetic field. However, for a self-gravitating gas the magnetic field 
is no longer "external" , but generated by the electron magnetic moments themselves, 
hence the relation between these variables cannot be controlled: it emerges from the 
dynamics of the gas and depends on the interplay of initial conditions. As we show in 




Figure 7. The figure displays the plot of magnetization M for initial 
temperatures </>(0) = 10 -3 , 10 -4 , 10~ 7 and Tm3i Tm4 (vertical axis in the left 
hand side) and Vmt (vertical axis in the right hand side) vs magnetic field 0. 
We considered initial magnetic field /3(0) = 5 X 10 — 5 and a positive value for the 
initial shear component along the z axis (the same value as in the plot of S + H 
in figure The plots show that larger values for the magnetization correspond 
to lower temperatures. The case with zero initial shear is displayed in figure [5] 



It is important to recall that the study we have undertaken here is based on 
an equation of state and a thermodynamical potential (f2) that come from a one- 
loop approximation. The behavior of the system in the regime (T 2 <C eB <C m^) 
may seem strange: the energy density does not increase with the temperature All 
this could follow as a consequence that the one-loop approximation may not be 
appropriate in some regimes for a magnetized system in presence of finite density 
and temperature. This issue deserves a separate study, possibly in the context of non- 
perturbative calculations at high magnetic fields, which explore a different approach 
that may be applicable to the magnetized gas that we have examined here (see [22] 
and references quoted therein). 

7. Conclusion 

We have examined the dynamical and thermodynamical behavior of a magnetized, 
self-gravitating electron gas at finite temperature, taken as the source of a simplified 
Bianchi I space time represented by a Kasner metric, which is the simplest geometry 
that allows us identify the magnetic field as the main source of anisotropy. We 
regard this configuration as a toy model that roughly approximates a grand canonical 
subsystem of a magnetized electron source in the conditions prevailing near of the 
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center and rotation axis of a compact object (in which spatial gradients of physical 
and kinematic variables may be regarded as negligible). The resulting Einstein- 
Maxwell field equations were transformed into a system of non-linear autonomous 
evolution equations, which were solved numerically in the collapsing regime for a 
chemical potential: /j, = 2m e , a magnetic field in the range B ~ 10 7 — 10 8 G and 
temperatures T ~ 10 3 - 10 7 K. 

For all initial conditions that we considered the gas evolves into a collapsing 
singularity, which can be (depending on the initial conditions) isotropic ( "point-like" ) 
or anisotropic ( "cigar-like" ) . We found that for lower initial values of the magnetic 
field and temperatures (B ^ 10 8 G, T ^ 10 3 K) the resulting singularity is always 
point- like singularity, independently of the initial values of the shear. This result may 
be connected with the stability of volume elements in less idealized configurations, 
an issue that is outside the scope of this paper and deserves a proper examination 
elsewhere. 

The collapse time decreases as the initial magnetic field increases, but we did not 
find a direct proportionality relation between this time and the initial temperature 
(this may be a consequence of having assumed a high density regime) . The behavior 
of the temperature as the collapse proceeds also depends on initial conditions: for 
initial temperatures T ~ 10 3 — 10 6 K the gas cools down and for values T > 10 7 K the 
temperature increases. 

Using the numerical solutions of the system we found a monotonic relation 
between the magnetization and the magnetic field and between the energy density 
and the temperature for high temperature values, but relation does not occur for low 
temperature values. This difference in behavior can be explained by the predominance 
of the magnetic effects at low temperatures, though it may also be due to the 
limitations introduced by assuming an equation of state and a thermodynamical 
potential based on a one-loop approximation. Looking at this issue is beyond the 
scope of this paper and will be pursued in a separate work. 

The study we have presented can be readily applied to examine hadronic systems 
(complying with suitable balance conditions and adequate chemical potentials). The 
methodology we have used can also serve as starting point to study the origin and 
the dynamics of primordial cosmological magnetic fields. These potential extensions 
of the present work are already under consideration for future articles. 
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Appendix A. Transformation of the integrals 

The definitions of T«, T±, Tu, F n are given by equations (j6|)-([8]), together with the 
coefficients ff9])- (| 12|) . Hence, their derivatives take the following form: 
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where the coefficients C5, Cfe, C7 are given by: 
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while h(x) = l/ (4cosh 2 (x/2)). The integrals (|5j)-([I2)>. (|A~7| - (jX9)) take the form: 
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Integrating by parts and making a change of variables v — + x 2 — /j]/<j> we finally 
obtain: 




v + /i) 2 — 1 ] h (u) dv , i = l,4, 

c > <"■/') = { )' < _ (A.H) 

v + /x) — 1 I H (u) dv , i = 5,7, 

with 

Fi(x)=fZf(v)dy, ft(i/)=[4coBh 2 (i//2)]-\ (A.12) 
References 

[1] Stuart L.Shapiro, Saul A. Teukolsky. 1983 "Black Holes, White Dwarfs, and Neutron Stars" ed 
John Wiley & Sons, Inc. 
S. Chakrabarty, Phys. Rev. D 43 (1991) 627. 
S. Chakrabarty, Phys. Rev. D 54 (1996) 1306 |hep-ph/9603406] . 

M. Chaichian, S. S. Masood, C. Montonen, A. Perez Martinez and H. Perez Rojas, Phys. Rev. 
Lett. 84 (2000) 5261 |hep-ph/9911218] . 

C. Y. Cardall, M. Prakash and J. M. Lattimer, Astrophys. J. 554 (2001) 322 [astro-ph/0011148] . 
R. G. Felipe, H. J. Mosquera Cuesta, A. Perez Martinez and H. Perez Rojas, Chin. J. Astron. 

Astrophys. 5 (2005) 399 |astro-ph/0207150] . 
A. Perez Martinez, H. Perez Rojas and H. Mosquera Cuesta, Int. J. Mod. Phys. D 17 (2008) 

2107 |arXiv: 071 1.09751 [astro-phi 1 . 
L. Paulucci, E. J. F errer, V. de la Incera and J. E. Horvath, Phys. Rev. D 83 (2011) 043009 

| arXiv:1010."304Tl [astro-ph.HE]]. 
A. P. Martinez, H. P. Rojas and H. J. Mosquera Cuesta, Eur. Ph ys. J. C 29, 111 ( 2003). 
R. G. Felipe and A. P. Martinez, J. Phys. G G 36 (2009) 075202 |arXiv:0812jj337l [astro-ph]]. 

E. J. Ferrer, V. de la Incera, J. P. Keith et at, Phys. Rev. C 82, 065802 (2010). 

R. G. Felipe, D. M. Paret and A. P. Martinez, Eur. Phys. J. A 47 (2011) 1 arXiv:1003."3254l 
[astro-ph.HE]]. 

A. Ulacia Rey, A. Perez Martinez and R. A. Sussman, Gen. Rel. Grav. 40 (2008) 1499 

| arXiv:0708.0593l [gr-qc]] . 
A. Ulacia Rey, A . Perez Martinez and R. A Sussman, Int. J. Mod. Phys. D 16 (2007) 481 
|gr-qc/0605054] . 

D. Manreza Paret, A. Perez Martinez, A. Ulacia Rey and R. A. Sussman, JCAP 1003 (2010) 
017 ; arXiv:0812!2508l [gr-ac]]. 

H. Y. Chiu, V. Canuto, and L. Fassio-Canuto, Phys. Rev. 176 (1968) 1438; V. Canuto and 
H. Y. Chiu,Phys. Rev. 173 (1968) 1229; V. Canuto and H. Y. Chiu, Phys. Rev. 173 (1968) 
1220; V. Canuto, and H. Y. Chiu, Phys. Rev. 173 (1968) 1210 
M. Chaichian, S. S. Masood, C. Montonen, A. Perez Martinez and H. Perez Rojas, Phys. Rev. 

Lett. 84 (2000) 5261 |hep-ph/9911218] . 
Charles. W. Misner, Kip. Thorne, John. Archibald. Wheeler, "Gravitation," Edit: 

W. H. Freeman NY, 1998. 
J. Wainwright, G. F. R. Ellis, "Dynamical system in cosmology," ed Cambridge University Press, 
1997. 

D. Ebert, K. G. Klim enko, M. A. Vdovichenko and A. S. Vshivtsev, Phys. Rev. D 61 (2000) 
025005 j hep-ph/9905253 ; D. Ebert and K. G. Klimenko, Nucl. Phys. A 728 (2003) 203 
|hep-ph/0305149| . 
C. O. Dib and O. Espinosa, Nucl. Phys. B 612 (2001) 492. 

A. Ayala, A. San chez, G. Piccinelli and S. Sahu, Phys. Rev. D 71 (2005) 023004 
|hep-ph/0412135]. 



